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1. INTRODUCTION 

Frequency deviation in Power System due to variation between generation and load shall be 
rectified within a fraction of seconds resulting in stability and security. Load Frequency Control (LFC) of an 
extensive power framework can be alluded as the issue of controlling the recurrence by directing the created 
units with reaction to change in stack [1]. For framework soundness, LFC must furnish recurrence with zero 
enduring state mistakes and tie-line trade varieties, high damping of recurrence motions and diminishing 
overshoot of the unsettling influence. The objectives specified are conveyed effectively in past works by 
various creators utilizing Fuzzy rationale PI and PID controllers [2, 3], ideal control [4, 5]. Variable structure 
control [6, 7], versatile and self-tuning control [8, 9]. Down the line, different tuning rules have picked up the 
consideration for the previously mentioned goals in which Internal Model Control (IMC) [10] is one among 
them. The LFC PID controller configuration utilizing Laurent arrangement is clarified by Padhan and 
Majhi [11]. Double PI controller tuning utilizing swam enhancement calculation is introduced in [12]. The 
two-degree-of-freedom internal model control scheme suggested by Tan [10] consists of two controllers with 
two tuning parameters where simultaneous tuning of the two parameters is difficult. In practice, a simple 
control structure with a fewer number of tuning parameters is desirable. The proposed control structure (see 
Figure 1) for LFC design consists of only one controller (Ge). Kasireddy et.al designed a PID controller for 
LFC through reduced model order [13]. 
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Figure 1. Proposed control structure 





In Figure 1, G and Gme®™ represent the power system dynamics and its model, respectively. For 
LFC, controller design is inconvenient because G results in higher order plant models, which are 
approximated by lower order transfer functions with time delay using a relay-based identification method. 

This paper has been alienated into 6 sections. Modeling of power system dynamics with necessary 
derivations discourse in section 2. In section 3, the PID controller design method is discussed followed by 
Section 4 in which the simulation results are presented. Section 5 deals with Robustness analysis and 
performance of a power system using Kharitonov’s rectangles followed by conclusions in section 6. 


2. MODELING OF POWER SYSTEM DYNAMICS 

Figure 2 shows single area power systems with a linear model. From Figure 2 it can be noticed that 
the power is supplied to the single area by a single generator. There are two types of turbine used for a 
generation: (a) non-reheated (NRT) and reheated (RT). 


Droop Characteristic 





Figure 2. Single area power system 


The plant model used for LFC without droop characteristics is 
G =G,G,G, (1) 


Where Gg, Gt, and Gp are the dynamics of the governor, turbine, and load & machine, respectively. For a 
reheated turbine, 


_ cIls+1 
‘ (T,s+1)(T,s +1) 


Where T, is a constant and c is the portion of the power generated by the reheat turbine in the total generated 
power. For non-reheated turbine T, = 0. The plant model used for LFC with droop characteristic is 


= G,G,G, 
G= ——— > _ (2) 
1+ (G,G,G, /R) 


From (1) and (2) can be represented by the second-order transfer function model 
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State space equations in the Jordan canonical form become 
x(t) = Ax(t)+ bu(t- 6, ) (4) 
y(t) = cx(t) (5) 
Where 
7 of op 
Acn sih ee [1 -1] 
0 = 1 =k 


When a relay test is performed with symmetrical relay of height th, then the expression for the limit cycle 
output for O < t < Om is 
y(t) =ce“x(0)+ cA (e* = I)bh (6) 
Let the half period of the limit cycle output be t. Then the expression for the limit cycle output for 


Om stSTis 
(7) 


y(t) = ce“ )x(9, ) -CA7 (e^ —I} bh 


The condition for a limit cycle output can be written as 
(8) 


y(0) = cx(0) = -y(t) = 0 


Substitution of t = T in (7) and use of (6) gives the initial value of the cycling states 


x(0) = (1+ e*) A7 (268%) -e^ —I)bh (9) 


When tp is the time instant at which the positive peak output occurs and tp 2 Om, then the expression of the 


peak output Ap becomes 
(10) 


Acs c(e^ )x(6, J-A (e^ _ I)bh) 


p 


and the expression for the peak time becomes 


T ee 
In 
T, = T, —tIT> 


(11) 





t,=9, + = 
1+e 


Substitution of A, b, and c in (9) and (10) give 


T, (1 pert )(2e% WT goth _ 1) _ T, (1 set )(2e (On )/T2 _ eT 1) -0 (12) 


wel 2 (13) 


A, =kh al +e Jie (14e UR ih — 4 
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The (11-13) are solved simultaneously to estimate Om, T1 and T2 from the measurements of t, Ap 
and tp. The relentless state gain k is thought to be known from the earlier or can be assessed from a stage flag 
test. Care has been taken to explain the arrangement of non-direct conditions, so intermingling may not occur 
to a false arrangement. 


3. PID CONTROLLER DESIGN 
The nominal complementary sensitivity function for load disturbance rejection can be obtained as 


1+ GG, 


To reject a step change in the load of the power system, the asymptotic constraint should be satisfied 
so that the closed loop internal stability can be achieved [3]. 


lim (1-T)=0 (15) 


The desired closed-loop complementary sensitivity function is proposed as 


(as? +0,8+1) 5 
LS = pr (16) 
(Bs+1) 


Where f is the only tuning parameter for obtaining the desired performance of the power system. As 


there always exists a trade-off between the nominal performance and robust performance, B must be tuned 
according to the desired choice. a; and a2 can be obtained from (15) and the constraint as 


T ( - P) eh j -T ( - P) eh 1 
T, T; 


T, =i, 


and 


4 
G27; E en’? 1/4 T, (17) 


2 


Using (14), (15) and second order Pade’ approximation for the time delay term, we get 


G, - 6(a,8” + 0,8 + 1)(1,87 +Ls +1) (18) 
ema Me a) 
m, m, 
Where 
(6T, +6T, +40, ) | (6T,T, + 4T,0,, + 4T,6,, + On ) 
1 = ———— 2 = a O UOUE 
6 6 


m, = 60, +248 - 6q, 
m, = 36ß° +0 — 6a, +1680. + 20,0, 
m, = 24870, + 20,0, +4B0* + 24p° 


The (18) can be written in the form of a PID controller with lead/lag filter as 
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2 
G, =K, gts en el (19) 
S 


: b,s°+b,s+1 





Where 
6a a 
K, = ' T, = Q, T, = 2 
km, ez 
m m 
a, =l, a, =l, b, = 2 b, = 
M, M, 


4. SIMULATION RESULTS 

Consider a power system with a non-reheated and a reheated turbine whose model parameters are 
given by Kp = 120, Tp = 20, Tr = 0.3, Tg = 0.08, R = 2.4, T; = 4.2 and c = 0.35 [11]. The identified models 
and controller settings (see Table 1) for the power system with non-reheated and reheated turbines are 
obtained using (11-13, 17). The Nyquist plots of the identified and actual models are shown in Figure 4 to 
illustrate the accuracy of the identification method. To get stable and robust response, B values in Table 1 are 
obtained from extensive simulation studies. Figure 3 and Figure 5 show the frequency change of the power 
system following a load demand AP4 = 0.01. The stability robustness is tested by changing the parameters of 
the system by 50%. From the simulation results, it is evident that the proposed method gives significantly 
improved performances than the Tan’s method. 


Table 1. Control parameters for identified model 


Model Type Identified Model Control Parameters 
NRT (WD) 12060" K.=2.0245, T=0.5005, Ti=0.1332, a;=29.0238, 
(28.4952s + 1)(0.2202s +1) ao= 15.1661, bj=28.6982, b2=5.77239, B=0.01 
NRT (D) 250e °s K.=0.7192, T;=0.2075, Tg=0.1159, a;=0.9212, 
2 028s? +12.765s +106.2 a2=0.1411, bj=0.1515, b2=0.0234, B=0.07 
RT CWD) 120e°°*" K.=3.6549, T=0.5797, Ti=0.2355, a\=24.4801, 
(23.2137s + 1)(0.9057s +1) ay=29.7725, b=24.0884, bo=20.2681, B=0.01 
E 235.3e 229s K.=1.0619, T=0.2107, Ty=0.1828, a;=1.154, 
(D) 1.79s? +16.9s +100 a:=0.1323, b\=0.1973, b2=0.0231, B=0.065 
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(b) Parameters of the system change by 50% 


Figure 3. Frequency deviation of the closed loop system with non-reheated turbine 
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Nyquist Diagram 
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Figure 4. Nyquist plots for the power system with non-reheated turbine 


(a) Nominal sytems 
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Figure 5. Frequency deviation of the closed loop system with reheated turbine 


5. ROBUSTNESS ANALYSIS AND PERFORMANCE 
In this section, Robustness of the system has been analyzed using Kharitonov’s Theorem. Closed- 
loop characteristic equation A,, (s)and denominator of the closed-loop transfer function T(s) are the 


polynomials that make the control system stable. Considering the forward-path and feedback-path transfer 
functions G(s) and H(s), characteristic equation is 


Aa (8) = 1+G(s)H(s) = 0 
Ag. (S)=a,S"° +a, 48" +...... + 84S +89 (20) 


For simplicity, assume that the leading coefficient an is constant and the coefficients have been normalized so 
that an = 1. The polynomial coefficients can then be expressed as 


aefarar™ |, i=0,1.......n-1 (21) 


(22) 
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According to Kharitonov’s Theorem, an nth-degree interval polynomial family described by 
(1a) and (1b) is robustly stable if and only if each of the four Kharitonov polynomials is stable, that is, all the 
roots of those polynomials have strictly negative real parts. 


For the system 


120 


O'S) = (6 0854 1)(0.3s +1)(20s +1) 


The characteristic equation is 
Ag (S) = 0.48s° + 7.624s* + 20.388 +121=0 (23) 
For +10% variations in the coefficients of the polynomial, The intervals of the polynomial will become 
a, € | 0.528,0.432 | a, € | 8.3864,6.8616 | a; € | 22.418,18.342 | a € [1 33.1,108.9 | 


Figure 6 shows Kharitonov’s rectangles rotate around the origin in a counter-clockwise direction to 
satisfy the monotonic phase increase property of Hurwitz polynomials. For clarity, the graph is zoomed in 
Figure 7 to show the zero-exclusion point. As the Kharitonov’s rectangles do not pass through the origin, it is 
concluded that the closed loop system guarantees the robust stability. 
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Figure 6. Response of the system with +10% variations in the coefficients of 
polynomial on the complex plane 
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Figure 7. Response of the system showing zero exclusion point 


Enhanced performance of PID load frequency controller for power systems (Dola Gobinda Padhan) 


124 o ISSN: 2252-8814 


Table 2 shows it is observed that the proposed method gives less Integral Absolute Error (IAE), 
Integral Squared Error (ISE) and Integral Time Absolute Error (ITAE) as compared to Tan’s method, so 
closed loop performance is improved. If we compare the total variations of the control signals, the results of 
both methods are almost same. Thus, with same control signals, the proposed method gives comparatively 
fewer errors. 


Table 2. Various errors and total variations 


Integral Time Integral Absolute Integral Squared a 
Type of Model Ne n c FGO & Pace 8 a Total variations 

NRT (WD)-Tan 0.5164 0.1061 0.001494 0.0199 
NRT (WD)-Proposed 0.5147 0.1008 0.001184 0.0127 
NRT (D)-Tan 0.5303 0.106 0.001363 0.0214 
NRT (D)-Proposed 0.521 0.1007 0.001175 0.0189 
RT (WD)-Tan 2.002 0.2007 0.002254 0.0735 
RT (WD)-Proposed 2.058 0.2061 0.002273 0.07906 
RT (D)-Tan 2.006 0.2007 0.01223 0.0729 

RT (D)-Proposed 1.966 0.2006 0.002264 0.0797 


6. CONCLUSION 

The Load Frequency Characteristics of a single-area power system with non-reheated and reheated 
turbines have been deliberated. The proposed method is flexible and gives satisfactory performance in 
nominal as well as the perturbed case. The proposed PID controller with a new control structure and a single 
tuning parameter ($) gave better performance than Tan’s controller. By showing the zero exclusion point by 
Kharitonov’s rectangles, it guarantees the robust stability for closed loop power systems. The proposed 
scheme can easily be extended to multi-area power systems. 
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